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We investigate numerically the clustering behavior of a system of phase oscillators with positive 
and negative couplings under a periodic external driving field with a bimodal distribution of driving 
phases. The phase distribution and the mean speed of the traveling state, as well as the order 
parameter for synchronization, are computed as the driving amplitude is varied. We observe that 
the periodically-driven system can also host traveling states for parameters in the same range as 
those for the case of a system without a driving field. The traveling speed is found to depend non- 
monotonically on the driving amplitude. In particular, oscillators divide into four clusters and move 
in pairs. Further, depending on the driving amplitude, two kinds of traveling mode arise: pairs of 
clusters traveling in the same direction (symmetric mode) and in opposite directions (antisymmetric 
mode). In the latter case (antisymmetric traveling mode), the average phase speed of the whole 
system apparently vanishes. A phenomenological argument for such behavior is given. 

PACS numbers: 05.45.Xt 


I. INTRODUCTION 

A system of sinusoidally-coupled phase oscillators, first 
introduced by Winfree pQ and later refined by Ku- 
ramoto [ 112 ], serves as a typical model for collective 
synchronization. When the coupling strength between 
oscillators is increased sufficiently, the system becomes 
ordered, forming a cluster of oscillators that move to¬ 
gether on a phase circle. Though there have appeared 
many variations of the Kuramoto model [4], the follow¬ 
ing extension of the model attracts interest with regard 
to the formation of clusters and their motion: When both 
repulsive and attractive couplings are appropriately in¬ 
troduced in the model [5], two clusters may be formed 
with the phase split by a difference of less than i r ra¬ 
dians. These two clusters, maintaining the phase differ¬ 
ence, travel on the phase circle; this is dubbed a traveling 
state. Recently, the dynamics of the traveling state has 
been studied in terms of the phase speed of the clus¬ 
ters [6]. According to that study, the positive (attrac¬ 
tively coupled) cluster chases after the negative (repul¬ 
sively coupled) cluster indefinitely. 

Other extensions of the model include a periodic driv¬ 
ing field, which may give rise to periodic synchroniza¬ 
tion: The order parameter may oscillate if the coupling 
is strong enough [7]. When the driving strength is large 
enough, a synchronized state is usually followed by a 
desynchronized state in a single period, and the dura- 
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tion of the former depends on the parameters of the sys¬ 
tem. Recently, we have shown that two clusters may 
be formed when a system of globally-coupled oscillators 
is driven periodically, which can explain the periodic 
synchronization-desynchronization behavior [8]. 

In this work, we consider a model that incorporates 
both features of the two types of extension: a system 
of positively and negatively coupled oscillators under a 
periodic sinusoidal driving field with a bimodal distribu¬ 
tion of driving phases. We observe, via numerical cal¬ 
culations, that four clusters emerge and travel in paired 
motion. In particular, two modes of traveling exist, de¬ 
pending on the phase of the driving field: The two pairs 
travel either in the same direction (symmetric mode) or 
in opposite directions (antisymmetric mode). In the lat¬ 
ter mode, albeit the parameters are the same, the system 
appears not to be in a traveling state. 

This paper consists of four sections: Section II de¬ 
scribes the oscillator model and its dynamics. In Section 
III, numerical results and the phenomenology of the two 
types of traveling mode are presented. Finally, a brief 
summary is given in Section IV. 


II. NUMERICAL CALCULATION 

We consider a system of N oscillators, the ith of which 
has an intrinsic frequency An oscillator is described 
by its phase fa and is coupled globally to the other oscil¬ 
lators. The dynamics of such a coupled oscillator system 
is governed by the set of Langevin equations of motion 
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(i = 1, N): 

& = Ui - ^2 sin (<j)i - (j)j) + Iq cos (Qt + cq), (1) 

i=i 

where the intrinsic frequencies { 07 } are assumed to be 
symmetrically distributed according to the Lorentzian 
distribution g(uS) = ( 7 / 7 r)(u ; 2 -fy 2 ) -1 . The second term 
on the right-hand side represents sinusoidal interactions 
between oscillators with strength FQ, which is positive 
or negative depending on whether or not the oscilla¬ 
tor follows the mean field A in Eq. ©> below. Specif¬ 
ically, it is taken from the distribution T(K) = (1 — 
p)S(K—K n ) + p5(K—K p ), where K n (< 0) and K p (> 0) 
are negative and positive coupling constants, with p de¬ 
noting the probability that an oscillator has a positive 
coupling constant. The final term describes a driving 
field with period r = 2i r/U, where the driving phases, 
the a^s, are distributed according to the bimodal distri¬ 
bution f(a) = (l/2)[5(a)+5(a— 7 r)], independently of the 
natural frequencies 07 . In other words, the driving term 
consists of Io cos Ut and Iq cosU(£ + t/ 2 ), time shifted by 
half the period r/2 = 7 r/fi, with equal probabilities. 

In order to measure the synchronization of the system, 
we introduce the complex order parameter 

1 N 

* ss jy xy*'-*= Ae", ( 2 ) 

3 = 1 

where the magnitude A measures the synchronization of 
the oscillators with the average phase 6. The order pa¬ 
rameter defined in Eq. ([2| allows us to reduce Eq. 0 to 
a single decoupled equation: 

(j)i = uji — KiAsm((f)i — 6) + Iq cos (Qt + 07). ( 3 ) 

In this work we resort to numerical methods to investi¬ 
gate the system. Using the second-order Runge-Kutta al¬ 
gorithm, we integrate Eq. (|3| with the time step 5t = 0.01 
for the system size N = 2000. Initially (t = 0), the <^’s 
are taken to be randomly distributed between 0 and 2tt 
for all i. We fix the period of the driving field to be 
r = 5.12, with the positive coupling constant K p = 1 
and p = 0.6. Most of the results in this work remain 
qualitatively unchanged even when we choose other val¬ 
ues of p for which a traveling state emerges. 

After the initial transient behavior, the system reaches 
stationarity and we obtain the time series of the order 
parameter as well as the time evolution of the phase 
distribution. Here, we examine the populations in a 
finite number of regions in phase space. Specifically, 
we divide one cycle of the phase into 72 different re¬ 
gions, the kth of which is defined by the phase interval 
(/c7T/36—7 t/ 72, fc7r/36+7r/72] (modulo 27r), and we obtain 
the number of oscillators belonging to the kth region 
(k = 0,..., 71). We also obtain the mean traveling speed 
(w) averaged over 30 ensembles, with each undergoing 
its own time evolution. Here, the traveling speed w is 
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FIG. 1: (color online) Magnitude A and phase in units of 7 r, 
0/ 7T (modulo 2) of the order parameter versus time t for a 
system of N = 2000 oscillators with Io — 1, p — 0.6, K n — 
—0.5, and 7 = 0.01. Also plotted is the cosine function of the 
driving field. 


obtained in the following way: At each time, we calcu¬ 
late the average value of the phase velocities Wi = <pi and 
take the time average over 10 periods of the driving field. 
Then, we take the absolute value to get w. 


III. RESULTS AND DISCUSSION 

Figure [l] present the time evolution of the magnitude A 
and the phase in units of 7 r, 6 /it (modulo 2) of the order 
parameter in a system of N = 2000 oscillators for Iq = 1, 
p = 0.6, K n = —0.5, and 7 = 0.01. Also plotted is the co¬ 
sine function of the driving field. Without a driving field, 
as described in Ref. 5, the traveling state emerges when 
the magnitude of the negative coupling constant is less 
than unity (\K n /K p \ < 1), and nearly half the oscilla¬ 
tors are positive (p « 0.5). For small driving amplitudes, 
for instance Jo = 1, the system exhibits periodic syn¬ 
chronization, as revealed by the behavior of A in Fig. 
What is different from the traveling state without a driv¬ 
ing field is the appearance of ripples both in A and in #, 
reflecting the periodic driving field in the system. Pre¬ 
sented in Fig. [2] is the distribution of oscillators, i.e., the 
number n & of oscillators belonging to the k -th region, at 
nine successive times in one period of the driving field 
from t (lowermost) to t + r (uppermost). Here, clearly, 
the two clusters, the larger one consisting of positive os¬ 
cillators and the smaller of negative ones, emerge twice 
in one period of the driving field. Subsequently, each 
cluster divides into two clusters in accordance with the 
driving phase a. In other words, four clusters charac¬ 
terized by the sign of coupling and by the driving phase 
(see below) emerge. After half the period of the driving 
field, two positive/negative clusters merge and generate 
a traveling peak. 

As the driving strength is increased, the cluster behav¬ 
ior of the system becomes more complex: Two clusters 
move in a pair in the same direction with a phase split 
while two pairs travel in the same direction or in opposite 
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FIG. 2: (color online) Number rik of oscillators belonging 
to region k with phase 4> satisfying /c7r/36— 7t/72 < 0 < 
/c7r/36+7r/72 (modulo 27r) at nine successive times during one 
period in the same system as in Fig. [l] From bottom to top, 
the nine data sets describe the distribution of oscillators rik 
at times t + (n/8)r for n — 0,1, ..., 8, respectively. The lines 
between data points are merely guides to the eye. For clarity, 
the curves have been successively shifted upward by 100. 


directions for some values of Jo, depending on the initial 
configuration. Thus, in some cases, the clusters do not 
appear to travel at all, as illustrated in Fig. [3J which 
presents the time evolution of A and 0/tt (modulo 2 ) in 
a system of N = 2000 oscillators with Iq = 1.6, p =» 0.6, 
K n = —0.5, and 7 = 0.01 for two different initial configu¬ 
rations of phases. Compared with the behavior in Fig. [lj 
the oscillatory behavior of A is more pronounced, and 
the phase changes more rapidly in Fig. [3^a). However, 
the mean traveling speed, which may also be estimated 
from the slope of a line joining the local maxima of 9 , is 
apparently smaller than that in Fig. [l] 

On the other hand, Fig. |3jb) reveals that the behav¬ 
ior of the same system can be significantly different for 
a different initial configuration: The order parameter A 
is not periodic in time and the mean traveling speed is 
apparently zero, as the phase 0 remains almost constant 
except for sporadic and abrupt changes of sign. To con¬ 
firm the apparently aperiodic behavior of A, we display 
in Fig. [4] the phase distribution rik corresponding to the 
data of Fig.[3jb). Manifestly, the phase distribution does 
not vary periodically in time: A single cluster at time t 
splits into two clusters, which do not merge into a sin¬ 
gle one, but remain as two at time t + r. Moreover, the 
two pairs of clusters travel in opposite directions. This is 
made clearer by probing the time evolution of the average 
phases of the aforementioned clusters, which are calcu¬ 
lated from the order parameters for all clusters in the 
system in Fig. [ 3 ] and are displayed in Fig. [5] The aver¬ 
age phase of each cluster is labeled according to the signs 
of the coupling (K = K p or K n ) and the driving phase 
(a = 0 or 7 r). While in Fig.[5|a) the two pairs of clusters 



FIG. 3: (color online) Magnitude A and phase 6 in units of 7 r 
(modulo 2) of the order parameter versus time t for a system 
of N = 2000 oscillators with / 0 — 1.6, p = 0.6, K n — —0.5, 
and 7 = 0.01. Also plotted is the cosine function of the driving 
field. Data sets in (a) and (b) have been obtained from the 
same system, but starting from different initial configurations. 



FIG. 4: (color online) Number rik of oscillators belonging to 
region k at nine successive times, as labeled, with the data 
obtained from the same system as that in Fig.[3|b). The lines 
between data points are merely guides to the eye. For clarity, 
the curves have been successively shifted upward by 100. 


consisting of positive and negative oscillators with the 
driving phase a = 0 and with the phase a = 7 r travel in 
the same direction, the two pairs of clusters travel in op¬ 
posite directions in Fig.[5|b). Note that positive and neg¬ 
ative oscillators that belong to the same pair (the same 
driving phase <a) move in the same direction. We may, 
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FIG. 5: (color online) Average phase 0 of each cluster ver¬ 
sus time £, as determined from the order parameters of each 
cluster in the systems corresponding to Figs. |3ja) and (b), 
respectively. The indices K and a labeling clusters are given 
in the legend. 

thus, conclude that two modes of traveling are possible in 
a driven oscillator system: the symmetric mode (one-way 
motion) and the antisymmetric mode (two-way motion). 
Regardless of the traveling mode, both pairs do travel, 
and the phase split is always less than i r. The average 
phase velocity, however, vanishes for the antisymmetric 
mode because two pairs travel in opposite directions at 
the same traveling speed. 

Whenever the antisymmetric mode comes into play, 
the mean traveling speed remains relatively small as 
shown in Fig.[ 6 j which presents the mean traveling speed 
(w) and the average order parameter (A) versus Iq in a 
system of N = 2000 oscillators with p = 0.6, K n = —0.5, 
and 7 = 0.01. Here, data only for the symmetric mode 
have been used in the calculation of ensemble averages. 
Also shown is (A) versus Iq for the system with only pos¬ 
itive couplings (p = 1). We observe that (A) has local 
minima and exhibits an oscillatory behavior. Note that 
the data shown in Fig. [3] describe the case Iq = 1.6, which 
corresponds to the first local minimum in Fig. [ 6 j 

In order to understand this behavior, we consider a 
system of positive oscillators, i.e., oscillators interacting 
via only positive couplings. Figure [7] shows how the order 
parameter A evolves in time t for four different values of 
Jo, two of which (J 0 = 3.5 and 5.4) correspond to the local 
minima of A in Fig. [ 6 ] while the other two (Io = 2.4 and 
4.4) correpond to the local maxima. For the latter two 
values of Jo, as indicated by the maxima of A, the two 
clusters overlap or are located rather close to each other. 
This occurs when the driving term is small, i.e., the co- 



FIG. 6 : (color online) Mean traveling speed (w) and average 
order parameter (A) versus driving amplitude Io for a sys¬ 
tem of N = 2000 oscillators with p — 0.6, K n — —0.5, and 
7 = 0.01. Also shown is (A) versus io for a system with¬ 
out negative coupling (p = 1). Error bars denote standard 
deviations, and lines are merely guides to the eye. 
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FIG. 7: (color online) Magnitude A of the order parameter 
versus time t for a system of N = 2000 positive oscillators 
(p = 1) with 7 = 0.01 and Io = 2.4, 3.5,4.4, and 5.4. The 
cosine function of the driving is also plotted. The second 
(Io = 3.5) and fourth (J 0 = 5.4) curves correspond to the 
values of io at which the order parameter reaches local minima 
in Fig. [ 6 ] 

sine function of the driving field takes a small value. The 
order parameter A reaches its maximum before the cosine 
function changes its sign. The once overlapping clusters 
become separated again as time goes on. When the di¬ 
rection of the driving field is reversed, the two clusters 
overlap again and then separate. On the other hand, for 
the former values of Io (= 3.5 and 5.4), A barely reaches 
a maximum as the sign of the cosine function changes; 
therefore, the values of (A) are smaller than those of the 
latter cases that exhibit double peaks in the time evolu¬ 
tion. 

We now return to systems consisting of both positive 
and negative oscillators. Suppose that pairs of clusters 
are established in the antisymmetric mode by the initial 
conditions or during the transient time. Two positive 
(or negative) clusters travel in opposite directions. If 
the driving amplitude Iq takes one of those values cor- 
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responding to the local maxima in Fig. [6j these clusters 
will have sufficient time to interact and to get to the sym¬ 
metric mode. Recall the presence of double peaks of A in 
the system of positive oscillators shown in Fig.[7| signify¬ 
ing the existence of two clusters for some time. Because 
the clusters are not much separated during this time, the 
interactions between clusters via coupling begins to have 
an effect, so that the motion of the pairs becomes sym¬ 
metric after several periods of the driving field. When J 0 
takes a value near one of the local minima in Fig. [6j on 
the other hand, this is less likely because the two pairs 
of clusters do not have enough time to evolve into the 
symmetric mode. As a result, the system may remain in 
the antisymmetric mode. The symmetric and the anti¬ 
symmetric modes are expected to appear (almost) pe¬ 
riodically with Jo because the oscillators tend to follow 
the driving field; if the system is in the (anti-)symmetric 
mode at some value of Jo, then the system is also in the 
(anti-)symmetric mode at the other value of Jo increased 
by A I (in our case A I ~ 2), except that the oscillators 
move longer on a phase circle. The oscillatory behavior 
persists for smaller values of p with lower values of (w) 
(data not shown). 

Finally, the aperiodicity appearing in Fig. [3|b) can be 
understood as follows: While clusters in the symmetric 
mode are rather far apart in phase from each other, in 
the antisymmetric mode, the distance between two clus¬ 
ters often becomes small, which makes the interactions 
between these clusters considerable. Such effects acceler¬ 


ate or decelerate the motions of clusters, resulting in the 
aperiodicity. 


IV. SUMMARY 

We have numerically studied systems of oscillators in¬ 
teracting via positive and negative couplings under a pe¬ 
riodic sinusoidal driving field with the driving phases fol¬ 
lowing a bimodal distribution. Observed is the forma¬ 
tion of two pairs of clusters, each consisting of positive 
and negative clusters according to the phase of the driv¬ 
ing field. Depending on the amplitude of the driving 
field, two pairs of clusters travel in the same direction or 
in opposite directions during each period of the driving 
field, which we call the symmetric and the antisymmetric 
modes of traveling, respectively. We have also found that 
the mean traveling speed of the symmetric mode depends 
non-monotonically on the driving amplitude and that the 
traveling speed is small in the antisymmetric mode. Dis¬ 
cussions on how the symmetric and the antisymmetric 
modes arise are given. 
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